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Abstract 

In this paper, we study the 3D axisymmetric Navier-Stokes Equations with swirl. 
We prove the global regularity of the 3D Navier-Stokes equations for a family of large 
anisotropic initial data. Moreover, we obtain a global bound of the solution in terms 
of its initial data in some LP norm. Our results also reveal some interesting dynamic 
growth behavior of the solution due to the interaction between the angular velocity 
and the angular vorticity fields. 

1 Introduction. 

Despite a great deal of effort by many mathematicians and physicists, the question of whether 
the solution of the 3D Navier-Stokes equations can develop a finite time singularity from 
a smooth initial condition with finite energy remains one of the most outstanding open 
problems |6j. A main difficulty in obtaining the global regularity of the 3D Navier-Stokes 
equations is due to the presence of the vortex stretching, which is absent for the 2D problem. 
So far, the global regularity of the 3D Navier-Stokes equations has been obtained only for 
initial data which are small in some scaling invariant norm [TH1 E3 EES]- But the analysis 
for small initial data does not generalize to the 3D Navier-Stokes with large data. A more 
refined analysis which takes into account the special nature of the nonlinearities and the 
anisotropic nature of the solution near the region of blow-up seems to be needed. 

In this paper, we study the global regularity of the axisymmetric Navier-Stokes equations 
with large initial data that have anisotropic scaling. Let u 9 and u e be the angular velocity 
and vorticity components of the 3D axisymmetric Navier-Stokes equations. We consider 
initial data for u e and ui e that have the following scaling property: 

u 6 {r,z,0) = -^U {er,z), u e (r,z,0) = ^W (er, z), (1) 
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where r = a/x 2 + y 2 , 5 and e are some small positive parameters, and the rescaled profiles 
Uo and Wq are bounded in L 2p and L 2q respectively for some p and q with p = 2q. We note 
that these initial data are not small. In fact, we have 

||Uo||L2(R2 x [ 0il ])||Vu || L 2( K 2 x [ 0)1 ]) = -J-^: > 1, 

for e small, where u is the initial velocity vector (we use bold letters to denote vector fields 
throughout this paper). Thus the classical regularity analysis for small initial data does not 
apply to these sets of anisotropic initial data. 

In this paper, we prove the global regularity of the 3D axisymmetric Navier-Stokes equa- 
tions for both initial data given by ([I]) and by exploring the anisotropic structure of 
the solution for e small. We also obtain a global bound on ||M 6, || i 2 P and ||ti; e '|| L 2< 7 in terms of 
their initial data. Note that by using the scaling invariance property of the Navier-Stokes 
equations, our global regularity result also applies to the following rescaled initial data 

u \r,z,0) = -^U (rA), u d (r,z,0) = -I^ (r, , (2) 

and 

u\r,z,0) = -U (^^), J f (r,z,0) = ±W (j =g ,Zy ( 3 ) 

For the rescaled initial data ([2]) and (j3J), our analysis suggests that ||M 9 || i 2 P may experience a 
rapid growth dynamically due to the contribution from Uq. On the other hand, the growth 
of ||a/|| L 2q is more moderate. See Remark 2 and (|32l-(!33l) in Section 3 for more discussions. 
This phenomenon is similar to the solution behavior of the ID model derived by Hou-Li in 
[9] where they show that u e is essentially bounded by its initial data while u e can experience 
large growth dynamically due to the contribution from {u d ) z and Uq. 

Note that the parameters e in the initial data ([I])-© and 5 in ([3]) measure the degree of 
anisotropy of the initial data. If 5 = 0, then the initial data ([3]) become isotropic, i.e. 

, N 1 TT fx y z\ 
Uo[x,y,z) = -U . 

e Ve e e/ 

Our analysis breaks down when there is no anisotropic scaling in the initial data, i.e. 5 = 0. 
Clearly, if the analysis could be extended to the case of 5 = 0, one would prove the global 
regularity of the 3D axisymmetric Navier-Stokes equations for general initial data by using 
the scaling invariance property of the Navier-Stokes equations. It is interesting to note that 
by using an anisotropic scaling of the initial data, we turn the global regularity of the 3D 
Navier-Stokes equations into a critical case of 5 = 0. 

We would like to emphasize that our global regularity results are obtained on a regular 
size domain, IR 2 x [0,1], for initial data (OQ). In this sense, our results are different from 
those global regularity results obtained for a thin domain, Q e = Q\ x [0, e] with Q\ being 
a bounded domain in R 2 . We remark that the global regularity of the 3D Navier-Stokes 
equations in a thin domain of the form Q e has been studied by Raugel and Sell in a series of 
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papers [23j [2U [25] . They prove the global regularity of the 3D Navier-Stokes equations under 
the assumption that || Vuollf^n \ < Coin-. This is an improvement over the classical global 
regularity result for small data, which requires || Vu ||^ 2 (Q e ) < C*e [23J. We may interpret 
our global regularity result with initial data (j2J) as a result on a generalized thin domain. 
Note that the initial data given by (j2J) satisfy the following bound: ||Vuo||| 2 (n £ ) = Coe~ 5+25 
(here 6 > can be made arbitrarily small), which is much larger than the corresponding 
bound C In - required by the global regularity analysis of Raugel and Sell in (23J, [24j [25] . 

This paper is motivated by the desire to understand how the local anisotropic structure 
of the solution near the region of a potential singularity may lead to the depletion of vortex 
stretching, thus preventing the formation of a finite time singularity To fully characterize 
the local solution structure of a potential singularity would be an extremely difficult task. 
Such study would shed useful light into our understanding of the dynamic depletion of 
vortex stretching and the global regularity of the Navier-Stokes equations. There has been 
some encouraging progress along this direction recently. In particular, Necas, Ruzicka and 
Sverak [20] and Tsai [28] have ruled out the possibility of an isotropic self-similar blow-up 
of the 3D Navier-Stokes equations. More recently, Chen-Strain-Tsai-Yau [5] prove that for 
axisymmetric 3D Navier-Stokes equations, if the velocity field u satisfies a scaling invariant 
blow-up rate 

\u(x,t)\<C*/y/r* + (T-t), 

where r = a/x 2 + y 2 , then such blow-up is not possible at t — T. This to some extent ex- 
cludes the locally self-similar isotropic blow-up of the axisymmetric Navier-Stokes equations. 
A related result with a stronger assumption can be found in [12] for the 3D Navier-Stokes 
equations. 

There have been also some numerical evidences on the anisotropic scaling of the solution 
near the region of a potential blow-up of the 3D Euler and Navier-Stokes equations, see e.g. 
[HI [H ED El [221 HSl HOI HI]. In particular, the recent work of Kerr [T5] and Hou-Li [TQl EE] 
gives some detailed description of the anisotropic scaling of the vorticity field near the inner 
core region of the maximum vorticity for two slightly perturbed antiparallel vortex tubes. 
The motivation for considering the rescaled initial data (T5]) and ([3]) is to understand whether 
such locally anisotropic blow-up is possible for the Navier-Stokes equations. 

We remark that Hou-Li have recently proved the global regularity of the 3D axisymmetric 
Navier-Stokes equations with some large anisotropic initial data [9J. The results obtained 
in this paper complement those presented in [9|. For more discussions on the axisymmetric 
solutions of the Navier-Stokes equations, we refer to [I] and [5] and the references cited 
there. The 2D Boussinesq equations are closely related to the 3D axisymmetric Navier- 
Stokes equations with swirl (away from the symmetry axis). Recently, Hou-Li [8] and Chae 
[3] have proved independently the global existence of the 2D viscous Boussinesq equations 
with viscosity entering only in the fluid equation, but the density equation remains inviscid. 

The remaining part of the paper is organized as follows. We give the formulation of the 
problem in Section 2. In Section 3, we prove a useful property of Riesz operators in R 4 x T 1 , 
where T 1 is the one-dimensional torus with periodicity 1. This property is needed to prove 
our global regularity analysis of the axisymmetric Navier-Stokes equations with anisotropic 
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initial data. Finally, we will present and prove the main result of this paper in Section 4. 



2 Formulation 

Consider the 3D axi-symmetric incompressible Navier-Stokes equations with swirl. 

u t + (u • V)u = -Vp + Au, 

V • u = 0, (4) 
u| i=0 = u (x), x = (x,y,z). 

Let 

e r =(-Ao), efl=M,-,0), e 2 = (0,0,1), 

be three unit vectors along the radial, the angular, and the z directions respectively, r = 
\Ac 2 + y 2 . We will decompose the velocity field as follows: 

u = v r (r, 2, i)e r + M (r, z, t)e e + t> z (r, z, t)e z . (5) 

In the above expression, u e is called the swirl component of the velocity field u. The vorticity 
field can be expressed similarly 

u = -(u e ) z (r, z, t)e r + c/(r, z, t)e e + -(ru e ) r (r, z, t)e z , (6) 

r 

where us 9 = v r z — v z r . 

One can derive evolution equations for u e and u e as follows (see e.g. [T51 14"]). 

u\ + v r u e r + v z u 6 z = (V 2 - 1) u e - -rfu\ (7) 



J t + v*J r + v*J z = ( V 2 - \ )u° + l - {(u e f) z + V, 

l ry A j ry » ' a* ry 
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V 2 - = (9) 

where is the angular component of the stream function, v r and v z can be expressed in 
terms of the angular stream function tp e as follows: 

v r = ~^f, v^-Ur^l (10) 
oz r or 

and V 2 is defined as 

V 2 = d 2 + -d r + «9 2 . (11) 
r 

Note that equations ([7j)-([9]) completely determine the evolution of the 3D axisymmetric 
Navier-Stokes equations once the initial condition is given. 



4 



By the well-known Caffarelli-Kohn-Nirenberg theory [2] , the singularity set of any suitable 
weak solution of the 3D Navier-Stokes equations has one-dimensional Hausdorff measure zero. 
Thus, in the case of axisymmetric 3D Navier-Stokes equations with swirl, if there is any 
singularity, it must be along the symmetry axis, i.e. the z-axis. Therefore, we should focus 
our effort to understand the possible singular behavior of the 3D Navier-Stokes equations 
near the symmetry axis at r = 0. 

As observed by Liu and Wang in [18], any smooth solution of the 3D axisymmetric 
Navier-Stokes equations must satisfy the following compatibility condition at r = 0: 

u%0, z, t) = c/(0, z, t) = VAO, z, t) = 0. (12) 

Thus we can rewrite u, u> and ip as follows: 

u 6 (r, z,t) = ru\(r, z,t), to e (r, z,t) = rui(r, z,t), ip 6 (r, z,t) = ripi(r, z,t). (13) 

In [9], Hou-Li have derived the following equivalent system for (ui, u>i, ipi): 

(Ui) t + t/'(Wl)r + V Z (Ui) z = 2^1 2 Mi + (u Uz + Ux„ H ") (14) 

r 

{uJl) t + V T (uJi) r + V Z (u 1 ) z = ((«i) 2 ) z + (uJlzz + Ulrr H ~) (15) 

r 

-(i>lzz + tPlrr + = Wl, (16) 



where v r and v z are defined as 

0{rWl) - - 1 9 (rVO- (17) 



dz r dr 

We note that the incompressibility condition implies that 

(rv r ) r + (rv z ) z = 0, (18) 

which can be verified directly from (jTT|) . 



3 A Useful Property of Riesz operators on R 4 x T 

Before we present our global regularity analysis, we would like to state and prove two 
technical lemmas regarding the property of a Riesz operator on M 4 x T\ where T 1 is the 
one-dimensional torus with periodicity 1. We first recall the following weighted Calderon- 
Zygmund inequality for a singular integral operator with a weight function which is in the 
yip class (see Stein pp. 194-217 for details). Let K be a Riesz operator in IR n and 
w{x) be a weight in the A p class (see page 194 of for definition). One can extend the 
Calderon-Zygmund inequality for the singular integral operator with the integral having 
weight function w(x). Specifically, for 1 < p < oo, there exists a uniform constant C such 
that 

/ \K * f\ p w(x)dx < C p [ \f\ p w(x)dx 1 (19) 
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for all / G L p (R n ). 

In this section, we prove a similar property for a Riesz operator on K 4 x T 1 . Let A = 
T^ =1 d 2 . + d 2 . Assume that f(r, z) belongs to the weighted L P (M. 4 x T 1 ) space with a weight 
function w(x) = \. Let u = u(r, z) , r = \J x\ + x\ + x\ + scf, be the solution of 

- A« = /, (20) 

with periodic boundary conditions along the z direction with period 1. We will prove the 
following lemma: 

Lemma 1. Assume that w(x) = w(r) belongs to the A p class in M. 5 . Then there exists a 
uniform positive constant C > such that 

/ | V 2 u\ p w{x)dxi ■ ■ ■ dx n _\dz < C / \f\ p w(x)dxi ■ ■ ■ dx n _\dz, (21) 

provided that f and u are both L p integrable with the weight w(x) for 1 < p < oo. 

Proof of Lemma 1. Let A > and denote the one-dimensional torus with periodicity v 
by iT 1 . For (x u ■ ■ ■ , x 4 , z) G M 4 x define 

'u\x x ,...^z)=u{\x x ,...,\x.M, 
f x (x 1} ■ • • ,x 4 , z) = X 2 f(Xxi, • ■ ■ , Xx 4 , \z). 



It is easy to see that 



Au x = f x , (xi,--- ,x 4 ,z) G M 4 x -T 1 . (22) 

A 



We claim that there exists a uniform constants C*i and Ci independent of A such that 

||(VV)w(a;)p|| LP(IR 4 x i T i ) < C x \f x w{x)~v\ LV( ^ijx ) + C 2 \\u x w(x)p\\ LP{R 4 x i T i ) . (23) 

To prove this, we introduce a smooth cutoff function <p(z) which satisfies < <p(z) < 1, 
(f)(z) = for z < — 1 or z > 2, and <p(z) — 1 for < z < 1. For A G (0,1], we define 
</> A (z) = </>(A,2). A simple computation gives 

- A(0V) = A / A - 2V0 A ■ Vm a - u x A<f) X , ( Xl , • • • , x 4 , ^) G 1R 5 . 

Since w(x) belongs to the A p class, the weighted Calderon-Zygmund inequality in M n (see 
( TT9]) ) implies that 

||(V 2 m a )w(x)p|| Lp(M 4 X i T i ) < ||V 2 (0 A it A )«;(x)p|| iP ( R 5) 

< C|| (</> A / A - 2V0 A ■ Vu A - u x A<p x )w(xf P || LP(R5) 

< C(II/ A w(^)^IIlp(k4 x i t i) + A 2 ||w A w(x)?|| iP(K 4 x i T i ) ) 

+ ^ll(VV)«;(x)^|| LP(K 4 x i T1) , (24) 
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where we have used the following estimate (note that <fi = <j)(z), w(x) = w(r) and u(x) 
u(r,z)): 

||V0 A -V M A Mx)^|| iP(R5) <CA( / \d z u x \ p w(x)dxy 



<C\[ I w(r)r 3 dr j \d z u x \ p dzy 



X 



< C p A 2 ||w A 'U;(x)p|| iP(R 4 x l T i ) + ^||(9V)^( a; )^lLP(M4xlTi)- 

In the last step of the above estimate, we have used the one- dimensional Sobolev inequality 
along the ^-direction: 

II?/ Il p < S^lli/P I 1 \\f) 2 v \\P 

for some positive constant C p depending on p. Thus, one has 

||(VV)u;(x)i|| iP(R 4 x i T i ) < C(||/ A -w(x)p HiP^s) + A 2 ||u A «;(a:)?|| L p ( || l 4 x i T i ) ) (25) 

for some uniform positive constant C > 0. For A > 1, we denote iV by the biggest integer 
less that A and compute 

- A(0m a ) = <pf x - 2V0 • Vw A - m a A0, (xi, • • • , x A , z) e M 5 . 

Similarly, by using the weighted Calderon-Zygmund inequality in M 5 , we obtain 

N\\(V 2 u x )w(x)p\\ LP{R , x i Tl) < C\\V 2 {(pu x )w{x)p\\ LPm 

< C\\(<Pf x - 2V0 ■ Vu x - n A A0) W (a;)F|| Lp{R5) 

< CN(\\f x w(x)p\\ LP{R4x i Tl) + \\u x w(x)p\\ LP{R4x i T1) ) 

+ yll(W)w(x)p|| LP(R 4 x i T i } , 

which gives 

||(VV)u;(x)p|| LP(R 4 x i T i ) < C(||/ A w(x)p || jL p C]R 4 x i T i) + ||^(x)*|| LJ , (R4x i T i)). ( 26 ) 



By ( l25l) and ( l26l) . we prove the claim ([23 
Note that 



\\(V 2 U X )w(x)p \\ LP{R 4 x l r i ) = A 2 v \\{y 2 u)w(x)v || iP(R 4 xT i } , 

||/ A w(x)p|| LP(R 4 x i T i ) = A 2_ p||/w(a;)p|| LP(R 4 xT i ) , 

1 5 1 

||m a w(x)p|| Lp(R 4 X i T i } = A"p||mw(x)p|| L p( R 4 xT i). 



Thus, we deduce from fl23l) that 



||(V 2 m)w(x)p|| lp(R 4 xT i ) < Ci||/w(x)p|| L p (R 4 )<t i ) + -t|||u«;(x)p||li»(r4 xT i). (27) 



By letting A — > oo in ([27)1 . we prove the lemma. 

Now we prove the following Lemma, which plays an important role in our global regularity 
analysis: 

Lemma 2. Assume that U\ and ipi are in L P (M. 2 x T 1 ) with 1 < p < oo, and ipi is the 
solution of (TT6]l with periodic boundary condition along the z-direction with period 1. Then 
we have 

IIV'I^IIlpQRSxT 1 ) < C||^l || Lp(R 2 xT 1 ) ) (28) 

for 1 < p < oo. 

Proof of Lemma 2. To prove ( 128|) . we observe that 

A _d^_ 3 _d_ 

dz 2 dr 2 dr 

is an axisymmetric Laplacian operator in five space dimensions, where r = \/ x 2 + x\ + x\ + x\ 
for x = (xi, x 2 , x 3 , x 4 , z) G R 5 . Let n = 5. Since —Aipi = u)\, Lemma 1 implies that 



/ 



\d^ipi\ p w(x)dx 1 ■ ■ ■ dx^dz < C \ \&i\ p w(x)dxi ■ ■ ■ dx^dz. (29) 

xT 1 JR^xT 1 

provided that w(x) = -j belongs to the A p class. We will show that w belongs to the A p 
class. In terms of the cylindrical coordinate, the above inequality fT2"9l can be written as 

1 poo \p / pi poo \p 

J \^ lzz \ p r 3 r- 2 drdz\ <C(j J \ui\ p r 3 r- 2 drdz\ , (30) 



which is exactly 

It remains to show that w(x) — \ belongs to the A p class. In the following, we will prove 
that w(x) = belongs to the A p class for < a < 4 (see page 194 of [26] for definition). 
To see this, let Br be a ball centered at x = with radius R. Denote by \Br\ the volume of 
Br, and p' be the conjugate of p, i.e. i + ^ = 1. We have 

w{x)dx ] [ — — - / w(x)~~dx 



\Br\ J Br J \\Br\ j Br 

n f pR pVW^ \ ( i pR pVW^ 1 

< ^R^ a (^R 5+a ^Y < en— = c. 



independent of R for < a < 4. Similarly we can prove the same estimate for Br centered 
at any point. This completes the proof of Lemma 2. 



4 The main result and its proof 

Now we state our main result in this paper. 

Theorem 1. Let < 5 < 1 be a given number which could be arbitrarily small. Assume that 

u e (r, z, 0) = -JzsU (er, z), u\r, z, 0) = ^W (er, z), (31) 

where U (r,z) = rUx{r,z), W {r,z) = rW x {r,~z), and U , Z7i E L 2p (R 2 x [0,1]), W , W x E 
L 2q (R 2 x [0, 1]) with p = 2q and q > 1/5. Further, we assume that the initial velocity field 
u E L 2 (R 2 x [0, 1]) and the initial conditions for ip e , u e and uo e are odd and periodic in 
z with period 1. Then there exists € (5,Uq,W ) > such that for all < e < e the 3D 
axisymmetric Navier-Stokes equations with initial data given by K31\) have a unique global 
regular solution. 

Remark 1. By using Theorem 1, we can easily obtain the global regularity of the initial 
data (El) by using a scaling argument. Specifically, if we denote by u the solution of the 
Navier-Stokes equations with initial data (EH), then w e (x, t) = \u ^) is the solution of 
the Navier-Stokes equations with the rescaled initial data (J2]). Thus the global regularity 
of the Navier-Stokes equations with initial data ([2]) follows immediately from Theorem 1. 
Using a slightly different rescaling, we can prove the global regularity of the Navier-Stokes 
equations with initial data ([3]). See Remark 4 at the end of this section for more detail. 

Remark 2. It is interesting to derive the corresponding global bound for the initial data 
([2]). Let {u e , u e ) be the solution corresponding to the initial data (|3Ti) . and {u e t1 uf) be the 
solution corresponding to the initial data (T5]). If we denote by u\ = u e e /r, to\ = to e /r, and 
ui = u e /r, uj\ = u e /r, then we have 

I r z t \ v z t 

u\(r, z, t) = -, -A u){{r, z, t) = 

Using the above scaling relationship and substituting / = u\ and g = uj\ into the global 
estimate (T4T1) . we obtain the following global estimate for u\ and uj\ with the initial data (TJ]): 

\\u\{t)\\L^ m < e-^+^WU.h^ + (2/^)^e- 2+5 / 2+1 ^||IU 1 ||i 29(CE) , (32) 
\H(t)\\L^ e ) < e-^ 2 «\\W x \\^ m + (C^e-^^mj^^, (33) 

where Q € = R 2 x [0, e]. Note that the first term on the right hand side of (|32l) is equal to 
|| u i(0)IU 2 p(n e ) and is smaller than the second term by a factor of 0(e s ^ 2 ). This implies that 
1 1 izj; 1 1 may experience a rapid growth dynamically for e small due to the contribution from 
the initial condition of u € . Similarly, the first term on the right hand side of fl33l) is equal 
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to ||^i(0)||i2 P (Q e -j, but is larger than the second term by a factor of 0(e~ 5 ). Thus the growth 
of \\UJ d e \\i2q is more moderate. This phenomenon is similar to the solution behavior of the 
ID model derived by Hou-Li in [5] where they show that (w l2 ) 2 + ujf satisfies a maximum 
principle. However, u\ can experience large growth dynamically due to the contribution from 
U\ z and u)\ at t = 0. 

Remark 3. It is worthwhile to point out that all the functions involved are in fact smooth 
for t > based on our assumptions and the dynamic control that we obtain. Here, the key 
is the dynamic control of ||wi(t) || L 2 P < C and ||wi(£)||lp < C. With these two estimates, the 
solution can be smoothly extended to all times. More specifically, the work of [2] and [17] 
implies that the velocity field u(x, t) is square integrable and smooth in the region where 
r > 1. Thus u(x, t) has a bounded L p norm in the region r > 1 for 2 < p < oo. On the other 
hand, we have \u e \ < C\ux\ and \w e \ < C\wi\ in the region where r is bounded. Therefore, 
our dynamic control of ||M 1 (t)||^2p < C (note that 2p > 3 by the assumption of Theorem 1) 
and < C lead to the following a priori bound on the gradient of u r e r + u z e z (see 

my- 

\\V(u r e r + M 2 e^)|| LP ( M 2 xT i ir<1 . 5 ) 

< C(||V x (u r e r + u z e z )\\ 

LP(K 2 xT 1 ,r<2) + ||V ■ (u T e r + U Z e z ) ||LP(K 2 xT 1 ,r<2) 

+ 1 1 u r e r + u z e z 1 1 L 2 ( R 2 x T i ir<2 ) ) 

< C(||a; 9 e e || L p( IR 2 xT i ir< 2) + ||M f e, r + n z e z || L 2( M 2 xT i ir< 2)) 

< oo, 

where we have used the fact that u> e eg = V x (u r e r + u z e z ) and V • (u r e r + u z e z ) = 0. The 
Sobolev embedding theory and the fact that u r e r + u z e z e L 2 (M 2 xT 1 ,r<2) imply that 
u r e r + u z e z G L 3 (IR 2 x T 1 ^ < 1.5) since p > |. Then the well-known regularity result (see 
[TB] ) implies that u(t) is smooth and LP integrable with p > 2 for t > 0. 

Proof of Theorem 1. By the assumptions of Theorem 1, we know that the initial velocity 
field u G L 2 (R 2 x [0,1]) and that u e \ t=0 G L 2p (R 2 x [0,1]), tu e \ t=0 G L 2q (R 2 x [0,1]) with 
p = 2q and q > 1/5. Using these properties of the initial data, one can show that there exists 
a finite time T > such that the axisymmetric Navier-Stokes has a unique regular solution 
for < t < T . In the rest of the proof, we will perform a dynamic estimate for the regular 
solution for < t < T and show that T can be made arbitrarily large. 

Denote Vt = R 2 x [0, 1]. Define / = |ux| p and g = \uJi\ 9 , with p = 2q. Multiplying (TJJJ by 
\ui\ 2p ~ 2 ui and integrating over Q, we get after using (fT8l) that 

/ f 2 rdrdz = / ip lz \u\\ 2p rdrdz + / \ui\ 2p ~ 2 ui(ui zz + U\„ H -)rdrdz. (34) 

2pdtJ n J n Jn r 
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Note that 



/ 

Jn 



ui\ 2p 2 uiU\ zz rdrdz 



-(2p-l) I \ Ul \ 2p - 2 (\ Ul \ z ) 2 rdrdz 
Jn 

-(2p-l) f (lu^- 1 ^,) 2 rdrdz 
Jn 



(2p-l) 

p2 



[ (f z y 

Jn 



rdrdz. 



Similarly, we obtain 



iti 



2p-2 \ ru lr)r , , 

p -Ux rdrdz 



\u\\ 2p 2 Ui(rui r ) r drdz 

\u 1 \ 2p - 2 {\u 1 \ r ) 2 rdrdz 



"(2p-l) / 
(2p-l 



(f r ) 2 rdrdz, 



n 



and 



Mi 



2p 2 Mi — —rdrdz 



\ui\ 2p 2 uiu lr drdz = - / (|wi| 2p ) r cirob 

= ~- f \ui\ 2p {0,z,t)dz = -- [ f 2 {0,z,t) 2 dz. 
V Jo V Jo 



Therefore, we have 

2p ~ 2 u\[u\ zz + U\ rr H —)rdrdz = 



ml 



3ui r N , , (2p — 1) 



1 r 

\Vf\ 2 rdrdz-- f 2 (0,z,t) 2 dz.(35) 
P Jo 



On the other hand, if we let p = 2q, q — p/{p — 1), we get 

ipufrdrdz < ||^iz||lp(Q)||/ 2 ||l«(Q) < llfelU^n) II/Il^) 



/ 



< C||a; 1 || L , ( n)||/ 2 |U,-(n) = Cy^U/ll*^ 



where we have used IIV'izIIlp(Q) < HV'izzll.LP(fi) which follows from the Poincare inequality along 
the z-direction, and the weighted Calderon-Zygmund estimate HV'i^IIlp(Q) < CII^iIIlp^)' 
which we have proved in Lemma 2. In order to apply Lemma 2, we need to show that 
■01 G L p (il). To see this, we use the fact that V r = — §-^ e G L m for any m > 2 and vp e is 
odd. Thus we can apply the Poincare inequality along ^-direction to show that i]r G L m (Q), 
which implies that ipi = ^ip G L S (Q) for any s > 1 in the region r > 1. In the region where 
r is bounded and < z < 2, the interior estimate of the general elliptic theory implies that 
V>i G L P (Q) whenever Ui G L P (Q). 
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Using the Sobolev interpolation inequality, 



L ^ <C\\f\\h\\^m a , (36) 



with a = 1 — ^ and q > 1, we obtain 



-- y Arfr^<C p ||^||2 2 ||/|| L2 2 1|V/||^ - i^-^IIV/lli,. (37) 



By the assumption, the initial conditions for ipi, U\, and U\ are odd functions in z. 
This oddness property of the solution is preserved dynamically. Thus we have /(r, 0, t) = 
g(r, 0,t) = 0. By using the Poincare inequality along the z-direction, we obtain 



\l 2 (u) < ||/z||z2(n)) IblU 2 ^) < lb«IU 2 (n)- (38) 
It follows from and ([37]) that 

^lX /2rdr ^- Cplbll ^ n)l|V/ " 2i2(n) " ^^ l|V/l1 ^) ' (39) 
Next, we multiply lu^l 29 ™ 2 ^! to ( fToT) and integrate over Q. We have 

*Q « r 9 

< / (P)M\ 2q -^irdrdz - ^^|| Vg\\l Hn) 
Jn Q 

1 /" 2 ! i (2o - 1) 

1 /2 

< (2 - -)|MUa ( n) ( / /^ 2(1 "^rdr^ - ^-=^|| V^||| 2(n) . 
Let p = p/2, q = -J^ = Then we obtain by using the Holder inequality that 

J fpg^ ^rdrdz^j < (^J f 2 rdrdz^j (^J g^^^p^rdrdzj 



< \\f\\h^ht r 



2 

iL 2 (n)iiy iiL 2 (n) ' 

where we have used the relationship p = 2q. Therefore, we obtain 



1 d f 1 2 i-- (2a — 1) 

2qdtj/ rdrdZ ~ ( 2 --)ll^ll^)ll/III 2 (n)ll^llL 2 (n)-^^l|V^I 



2 

L 2 (Q) 



4 



,2(i-|) (2g-l) , IT7 ,„ a 
2g 2 



< C!l/ll£ 2 ( n)!ML 2 (n) -^-^HV^L 2 (n) 



q 2p 2 ll/lli2 (^) + 4 g 2 W9\\Li(n) 2q2 HV0|| £ 2 (n) 
S O f ^i||v/||i, w) -^l)||V,||i, (m , (40) 
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where we have used the Poincare inequality (|38|) in the last inequality. 

Multiplying (|39|) by C q and adding the resulting equation to (|40|) . we get 



— ( C q [ f 2 rdrdz H / g 2 rdrdz] 

< c q (c^s||f 2(n) ||V/||| 2(n) - ^^l|v/||i 2(n) 



+c q ^^\\vf\\l 2m -^^\\v 9 \^ m 

<cJcM\im-^^)\\^f\\Un)- HI) 



^nyiii2 (n) - - ^ 2 - j ||v/||| 2(n) . 

Using ( BID , we have 

m(r,z,0) = e 5 f/ 1 (er,z), Wl (r ) z,0) = £ f lVi(er,z) ) (42) 
where C/i(f, 5) and Wi(f, 5) have bounded L 2p and L 2<? norms independent of e. We obtain 

IMlf 2 = e*~«||Wi||i*(n), (43) 
H/olli = e 5 ^||f/i||^(c). (44) 

We would like to choose e small enough to ensure that the right hand side of (|41|) is 
negative. To this end, we require that 



c P \\g\\h ia) (t)<^^-- («) 



If ( H5j) holds, we would get 



— — / f rdrdz + — / o 2 rdrd,z < 0. (46) 
This implies that 

— - / f 2 rdrdz H / g 2 rdrdz < — f f^rdrdz H / g 2 rdrdz. (47) 

2p 2g J n " 2p J n 2g J 

This gives a global bound on || Million) and ||o>i||.L2g(n) in terms of their initial data. In 
particular, we have by using (|43!) - (l44l) that 

\\9\\l m (t) < lbo|li ( a) + (C 3 /2)^||/ |l! 2( n) 

< e^WWxh^ + {C q /2)he 25 -\\\U l f L , P{sl) . (48) 



13 



Using (PfBl) . we conclude that the condition f|4"5]) is satisfied if we choose e to satisfy 

C P (^--llWill^n) + (C g /2)^6 25 -i||C/ 1 ||i 2p(n) ) < (49) 

which is true by taking q > 1/5 and e small enough. This condition would guarantee that the 
global estimate (14TI) is valid for all time. The global estimate (T4T1) gives a global bound on 
||'"i|U 2 p(r2) and ll^xll^^). Using this global estimate, we can easily obtain a global bound on 
||M 6 '||i2p and 1 1 uj 9 1 1 i,2q over a bounded domain with r < R, which proves the global regularity 
of the 3D Navier-Stokes equations (see Remark 3). This completes the proof of Theorem 1. 

Remark 4. The global regularity of the Navier-Stokes equations with initial data ([3]) follows 
almost exactly as the proof of Theorem 1 except that the bounds for ||<?o|U 2 and ||/o|U 2 are 
different. By using the same rescaling, we just need to consider the rescaled initial data of 
the form 

u e (r, z, 0) = U (eV, z) , J (r, z, 0) = W (e 5 r, z). (50) 
Due to the difference in the scaling, the bounds for H^olU 2 ? and ||/o||l 2< ? are now given by 

IMlt = e^-^llWiU^n), (51) 
WhWh = e 5(1 - 1/?3) ||£/i|U 2 ^). (52) 
Thus the condition (T4"5]) is satisfied if we choose q > 1 and e small enough to satisfy 

C P (e^WW^w + {Cj2)h^\\U4l^ < ^-r^, (53) 
Condition (153]) will guarantee that the global estimate ff4T|) is valid for all time. 
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